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ABSTRACT 

We present calculations of thermal evolution of Hot Jupiters with various masses and effective 
temperatures under Ohmic dissipation. The resulting evolutionary sequences show a clear tendency 
towards inflated radii for effective temperatures that give rise to significant ionization of alkali metals 
in the atmosphere, compatible with the trend of the data. The degree of inflation shows that Ohmic 
dissipation, along with the likely variability in heavy element content can account for all of the 
currently detected radius anomalies. Furthermore, we find that in absence of a massive core, low- 
mass hot Jupiters can over-flow their Roche-lobes and evaporate on Gyr time-scales, possibly leaving 
behind small rocky cores. 



1. INTRODUCTION 

Over the last decade, novel discoveries of transiting 
extra-solar planets have often been fraught by unex- 
plained radius anomalies. In particular, many gas gi- 
ant planets, residing on orbits in extreme proximity to 
their host stars, have been found to have much larger 
radii than what was previously thought possible in the 
context of "s tandard" gas giant theory (Stevenson"1982^ 
Guillot|1999 ). A number of possible explanations to this 
pro blein have been proposed , most notably tidal heat- 
ing ( Bodenheimer et al. [200 1', '2003), kinetic heating du e 
to breaking gravit y waves ( Guillot fc Showman 2002), 
enhanced opacity (Burrows et al.||2007 ), semi-convection 
( [Chabrier fc Baraf!e|20U7[ ) and turb ulent burial of heat 



by a mechanical greenhouse effect ( [Voudin fc Mitchell 
20101. However, it appears unlikely that any of the 



above solutions can simultaneously explain all o f the 
observed anomalies (Fortney fc Nettelmann 2010). Re- 
cently, a new Ohmic heating mechanism, that relies on 
the electro-magnetic interactions between atmospheric 
flows and the planetary magnetic field, has been sug- 
geste d as a promising explanatio n to the inflation prob- 
lem ( [Batygin fc Stevenson||2010[ ). 

A characteristic that is unique to hot Jupiters is the 
enormous amount of incident energy that they receive 
from their host stars. The extreme irradiation naturally 
results in high at mospheric tempera tures (sometimes in 
excess of 2000K) (Spiegel et al. 20091, and the latitudinal 
dependence of this heating leads to fast, super-rotating 
equatorial j ets with characteristic wind-speeds of order 
^ 1 km/s ( [Sho wman et al.||2009|). Consequently, the 
atmospheres of hot J upitcrs arc hot enough to thermally 
ionize alkali metals, that are present in trace abundances, 
and give rise to a finite electrical conductivity of order 
^ 10~^ — 1 S/m, depending primarily on the effective 
temperature. Rapid advection of the ions by the jets, 
while immersed in the planetary magnetic field, leads 
to an induction of an electro-motive force that sets up 
electrical current loops through the deep interior and the 
atmosphere of the planet. These electrical currents give 
rise to Ohmic heating as they flow throughout the planet. 

|kbatygin@gps.caltech.edu| 



It has been shown, in the context of a static structural 
model, that the level of resulting dissipation in the inte- 
rior is approxim ately that required to mainta in the radii 
of hot Jupiters (Batygin & Stevenson 2010). However, 
time-evolution ot Ohmicaliy heated gas-giants remains 
an open question, since the degree of dissipation exhibits 
strong dependence on the structure of the planet, and 
particularly on the planet's Teff^ mass, and the location 
of the radiative/convective boundary. More explicitly, 
the thermal structure can be affected by the heating and 
its radial variation as well as the elapsed time. At the 
same time, the conductivity arises from thermal ioniza- 
tion, causing the interior structure to feedback on the 
heating. Thus, the structural profile of the planet dic- 
tates the Ohmic heating distribution. As a result, self- 
consistent evolutionary calculations are required to as- 
sess whether Ohmic dissipation can indeed resolve the 
inflation problem. To perform such calculations is the 
primary purpose of this study. 

Here, we construct an approximate model with specific 
assumptions, that couples calculations of Ohmic heating 
and structural evolution, and show how this heating to- 
gether with the likely variability of heavy element abun- 
dances can provide an explanation for all of the currently 
observed radius anomalies. We compute the maximum 
radii attainable by evolved Hot Jupiters as functions of 
planetary mass and effective temperature and find a ten- 
dency for larger radii for effective temperatures that cor- 
respond to the onset of significant conductivity arising 
from alkali metals in the atmosphere. This behavior is 
compatible with the trend of the data. Furthermore, our 
calculations suggest that in absence of a relatively mas- 
sive, high-metallicity core, low-mass hot Jupiters can ex- 
pand beyond their Roche-lobes and spill their envelopes 
onto their parent stars on billion year time-scales, possi- 
bly leaving behind small rocky cores. 

The paper is structured as follows. In section 2, we con- 
sider the energetics of the Ohmic mechanism, and show 
that in steady state, the degree of dissipation is limited 
by the thermodynamic efficiency of the atmosphere. In 
section 3, we present a simple model for the damping 
of the global circulation by the Lorentz force, and show 
that characteristic efficiency of the Ohmic mechanism is 
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of order a few percent. In section 4, we describe the 
coupled Ohmic heating/themal evolution model that we 
employ in our calculations. In section 5, we present the 
theoretical curves that delineate mass-radius-Te// space 
and compare our results with observational data. We 
conclude and discuss avenues for future improvement of 
the Ohmic inflation model in section 6. 

2. WORK - OHMIC DISSIPATION THEOREM 

Let us begin by considering the energetics of the Ohmic 
dissipation mechanism. We start by writing out the invi- 
cid Navier-Stokes equation, where turbulent stresses have 
been neglected and the last term is the Lorentz force: 
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In the above equation, v is the velocity, fl is the rota- 
tional velocity, P is pressure, p is density, g is the gravita- 
tional acceleration, J is the current density and B is the 
background planetary magnetic field. There exists a vast 
literature on the subject of energy budgets of the vari- 
ous terms in t he N avier-Stokes equa t ion (s ee for example 
HoltonI ( [T992I ) and lPeixoto k OortI ( [1992] )). Although a 
discussion ot the global picture is important for under- 
standing the details of atmospheric dynamics, the aim of 
our study is limited to the Ohmic mechanism, so we shall 
focus our discussion on the Lorentz term. The power, or 
rate of change of kinetic energy energy of the fluid per 
unit volume, provided solely by the Lorentz force is 
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(2) 



If the quantity on the right is positive, the Lorentz force 
adds kinetic energy to the flow, while if it is negative, 
kinetic energy is drained from the flow. Recall that MHD 
Ohm's law reads 



/ = (t(vx B - V$) 



(3) 



where a is the electrical conductivity and $ is the electric 
scalar potential. Using a standard vector-calculus iden- 
tity, we can write v-JxB^—J-vxB. Consequently, 
the energy equation can be re-written as 



V- J X B = J-V$ 

(T 



(4) 



The first term on the RHS, is the Ohmic dissipation. 
Upon integration, by Gauss's theorem, the last term van- 
ishes, since we require no radial current at the outer 
boundary: 



J-S/<i>dV = 



V-(J$)dV^ = y)(J$)-da = 

As a result, we discover that in steady-state, Ohmic dis- 
sipation is work done by the flow: 
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(6) 



It is important to note, however, that the work done 
by the flow is limited by the efhciency factor i.e. the 
fraction of insolation that is available to do useful work. 



In practice, this means that the total Ohmic dissipation 
rate should be rather insensitive to the magnetic field 
strength, once the field is larger than some critical "sat- 
uration" value. Particularly, in the saturated case, the 
flow velocity should scale inversely with the magnetic 
field. Thus, in such a regime, changing the field will not 
change dissipation or the planet radius. Rather the effi- 
ciency factor plays the governing role. Although numeri- 
cal simulations are required to quantitatively understand 
the saturation field strength accurately, based on dimen- 
sionless number analysis, it is likely that the critical field 
strength is not overwhelmingly high (^ 1 Gauss or so). 

3. MAGNETIC DAMPING OF THE GLOBAL CIRCULATION 
AND THE EFFICIENCY OF THE OHMIC MECHANISM 

In the limit where the Lorentz force is dynamically neg- 
ligible, the wind speeds may be interpreted from global 
circulation models (GCMs). In this case, the computa- 
tion of Ohmic dissipation is straight-forward, although 
it is noteworthy that G CM's carry an intrinsic range 



of numerical uncertainty ( Cho et al. 2008 Showman et 



all 2010 Heng et al. 2011). I'or sufiiciently high atmo- 
spheric ionization levels and/or planets with strong mag- 
netic fields, however, the onset of flux-freezing can act to 
slow down the winds, effectively capping the degree of 
Ohmic dissipation. The total Ohmic dissipation can be 
written as some fraction of the absorbed sunlight; we 
refer to this as the efficiency, e. This efhciency must 
obviously go to zero as one goes to very low conduc- 
tivities or low background field strength and is limited 
from above by thermodynamic (Carnot efhciency) con- 
siderations. It is further noteworthy that only a part of 
the total Ohmic dissipation (-- 5% for Teff > 1400K; 
- 10-''% for Teff = lOOOK) is deposited in the region 
of relevance to our inflation mechanism (the convective 
interior) ; the rest is dissipated at higher level^ (see Fig- 
ure 4 and the associated discussion). The primary aim 
of this section is to derive a quantitative estimate for the 
Ohmic efficiency. 

The limitation to efficiency can arise in two ways: 
through the reduction of the wind speeds as a direct con- 
sequence of the Lorentz force or through the reduction 
of the equator to pole temperature difference because of 
meridional circulation. We have not performed full MHD 
global circulation mod els. Instead , following [Schneider] 
'fc Li ndzen] ( [1977] ) and [Held fc Hou| (|1980|, we develop 
a simple, analytical model tor the mid-latitude circu- 
lation in hot Jupiters, crudely accounting for magneto- 
hydrodynamic interactions. 

Let us begin by simplifying the Lorentz force. Since the 
electric field is of the same order of magnitude as w x B^i-p 

as defined by the boundary condition, we write J a{v x 
Bdip)- Collecting the conductivity, magnetic field and 
fluid density into a damping timescale and retaining only 
the component that opposes the flow, we have 
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In other words, the effect of the Lorentz force is to in- 
troduce a term into the equation of motion that can 

^ This is because for lower T^t f the atmosphere is characterized 
by higher resistance, and thus clissipates fractionally more energy. 
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be approximated as a Rayleigh drag with characteristic 
timescale 
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A damping time-scale of t he same functional form was 
used in GCM simulations of Perna et al. ( 2010 1 to mimic 
the the Lorentz force, in order to obtain an estimate of 
Ohmic dissipation in the atmosphere. As will be revealed 
below, this term introduces both direct drag as well as 
meridional circulation into an otherwise geostrophic so- 
lution. 

Next we assume steady state and azimuthal symme- 
try, so all time and zonal derivatives in the Navier-Stokes 
equation vanish. Finally, under an assumption of a con- 
stant Coriolis parameter (f-plane approximation), the 
horizontal Boussinesq equations of motion reduce to an 
Ekman balance: 

fVy ^ ^ (9) 



, 1 dP' Vy 

fvx = ^ 

Po oy Tl 



(10) 



where we have expressed all quantities in a local carte- 
sian coordinate system i.e. ^ v^-, vg ^ Vy, and P' is 
the component of pressure that is not compensated by 
gravity. Under the same assumptions, the vertical mo- 
mentum equation, neglecting Coriolis effects simplifies to 
the hydrostatic equation 
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where a = 1/ (0) is the coefficient of thermal expansion 
and e' is the potential temperature departure from the 
background state 6(z). Stable stratification is implicit 
in the problem and we write the heat equation in accord 
with Newtonian cooling: 
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(12) 



where 8'"'"' is the radiative deviation from the back- 
ground state, and tjv = CpP/AgasBT^jj is the Newto- 
nian cooling timescale. The physical meaning of 8''°'^ is 
that if no meridional flow is present to transport heat 
from the equator to the pole, the potential temperature 
profile takes the form 8 = 8-^8' = 8-^ 8™'^. 

To obtain a solution for the meridional circulation, 
let us introduce a stream-function ^(y, z) such that 



d'^ /dz and Vz 



-d'^/dy. Upon differentiating 



the hydrostatic equation with respect to z, differentiat- 
ing the y-momentum equation with respect to y and sub- 
tracting the first from the latter, we obtain a modified 
thermal wind equation. 
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Multiplying through by de/dz, expressing the velocities 
in terms of the stream-function, and applying the heat 
equation, we arrive at a differential equation for the po- 



tential temperature perturbation: 
^ fN^\ 1 d^e' 



2 oracZ 



d^e 



dz^ 



(14) 



where N'^ = g (8) ^ de/dz is the square of the Brunt- 
Vaisala frequency. Note that the Lorentz force is not 
necessarily anti-parallel to the zonal velocity so the core- 
spodence we use here is designed only to find the effect 
of the drag. Next we take the solutions to be of the 
functional form: 

^i^^a sin{T:y/R) sin(7rz/i/) 
Vx = -Vxo sin(7r?//i?) cos(7rz/ H) 
Qvad ^ Qvad cos{ny/R) sininz/H) 

e' ^e'oCOs{TTy/R)sm{7rz/H) (15) 

where H is the atmospheric scale- height, R is the plan- 
etary radius, and all terms with the subscript "0" are 
undetermined constants. Upon substitution of these so- 
lutions into the above equation, we obtain an algebraic 
expression for the actual potential temperature deviation 
as a function of the radiative equilibrium deviation and 
other system parameters. This expression can be sub- 
stituted back into the thermal wind equation to yield 
an order of magnitude equation for zonal wind velocity 
reduction due to the drag: 
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(16) 

where A8''°'' is the equator-to-pole potential tempera- 
ture difference that the planet would have if no merid- 
ional circulation was present. Consider the limiting cases 
of this equation. If we take the drag to be negligible 
(tl — J> oo), we have 8' = 8'""'', so the global poten- 
tial temperature becomes 8 = 6 -|- 8''"'', as Vy 0. 
Thus the wind velocity equation reduces to the standard 
thermal wind equation. If we consider the drag to be 
the dominant effect {tl 0), then the potential tem- 
perature perturbation again becomes equivalent to the 
radiative potential temperature perturbation, since all 
circulation effectively stops. This is a regime where the 
above treatment of the Lorentz force is highly inaccurate. 
The transition between the two regimes takes place when 

fTL 1. ^ 

The functional form for the stream-function that we 
introduced above satisfies the no-normal flow boundary 
conditions at the edges of the modeled "cell" . The cir- 
culation it implies, however, is somewhat unrealistic and 
is not globally observed on either terrestrial or solar sys- 
tem gas giant atmospheres. This is largely because in 
order to obtain the analytical solution presented above, 
we made a series of simplifying approximations and im- 
plicitly omitted any discussion of turbulent Reynolds 
stresses, V • (v'^v'y). In a detailed model for atmospheric 
circulation, the resulting eddy momentum fluxes are es- 
sential to clos ing the angular mome ntum cycle on gas 
giant planets ( Liu &: Schneider]|2010 ), although formally 
the above solution also closes Th e angular momentum cy- 
cle ( Schneider fc Lindzen|1977 |. Consequently, we stress 
that the above picture is not intended to be a quantita- 
tively good representation of the atmospheric dynamics 
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Fig. 1. — A comparison between interior current geometries, induced by (A) a single jet (B) two counter-rotating jets, each in one 
hemisphere and (C) a triple jet with retrograde equatorial flow. The small arrows inside the planet show the induced current (actual model 
output). Additionally, dipole magnetic field lines as well as graphical representations of the currents are presented. Atmospheric flows are 
represented as into-the-board and out-of-the-board arrows on the sides of the planets. 



on hot Jupiters. However, we do suggest that the model 
yields the correct behavior of the wind velocity reduction, 
to an order of magnitude. 

Direct application of the above formalism to our mod- 
els is not possible because of the large variation in con- 
ductivity and density. However, we can still use the re- 
sulting quantitative estimate of Ohmic efficiency as a 
guide, since any reasonable choices of parameters (e.g. 
the scale factors in equation 8), yield a Lorentz damp- 
ing timescale, r^, that is far longer than the rotational 
timescale, implying a limited effect on the zonal winds. 
In other words, for most hot Jupiter regimes, it is reason- 
able to approximate the effect of the Lorentz force as a 
Rayleigh drag. The efficiency of conversion of zonal wind 
energy into heating is then given by the ratio of the ki- 
netic energy dissipated by the drag term in the equations 
of motion to the incident stellar flux: 



er^ pv^H/{TLasBT'c 
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A value of this order is in accord with the results of 
simplified GCM simulations of ho t Jupiters w here mag- 
netic effects are modeled as drag (Perna et al. 2010). In 
order to qualitatively explore the ettect of variable effi- 
ciency, we compute three sets of models, characterized 
by 1%, 3%, and 5% Ohmic efficiencies. Note that at 
low {Teff ~ lOOOX) temperatures, such efficiencies are 
likely to be overestimates, because of very low conduc- 
tivity. However, the mechanism at very low conductivity 
is inactive anyway, so we retain the same efficiencies for 
consistency. From a computational point of view, a con- 
stant efficiency means that although the dissipation is a 
function of the radial coordinate, it is scaled such that 
the integrated Ohmic heating amounts to 1%, 3%, or 5% 
of the insolation. In other words, we only adjust the 
magnitude of the heating, not its radial distibution. 

4. COUPLED OHMIC HEATING/STRUCTURAL 
EVOLUTION MODEL 

In this study, the calculation of the Ohmic h eating 
was performed in effec tively the same way as in (Baty- 
gin & Stevenson 2010), but with the atmospheric tem- 



perature structure computed rather than assumed. In 
other words, as the planet evolves along its evolution- 
ary sequence, the conductivity profile everywhere in the 
planet is recalculated at every time-step. Naturally, the 
heating profile is updated at every time-step as well. A 
number of studies in the past have concentrated a ll of 
the (tidal) heating in the deep interior of the p lanet (Bo- 
denheimer et al.||2001[ [2003{ (ibgui et al.||201C)[ ). It is not 
known whether the tidal heating is deep-seated; it may 
be in the atmosphere and thus unavailable for inflation. 
By contrast the Ohmic model explicitly computes the 
relative amounts in the interior and atmosphere, so the 
actual calculated heating distribution is implemented in 
the structural calculations. Note that in our model, heat- 
ing of the adiabat ic interior as well as heating of the deep 
atmosphere (as in Guillot & Showman (2002)) contribute 
to the radial evolution. 

We solve the simplified steady- state induction equation 
(see Batygin & Stevenson (2010) for derivation) 



V • 0-V$ = V • cr (^t/ X BdipJ 



(18) 



All of these planets can be expected to have internal 
dynamos with field strengths similar to Jupiter {\B\ ^ 
few Gauss) or perhaps larger because of the lower den- 
sity in the dynamo generating region. It is implicitly 
assumed the magnetic field remains roughly dipolar at 
the planetary surface, in spite of the induced current. 
This assumption may be generally valid in most hot 
Jupiter atmospheres, as the magnetic Reynold's num- 
ber. Rem — VL/t], where rj is the magnetic diffusivity, is 
typically of order unity at most, rendering dynamo ac- 
tion unlikely in the atmosphere. However, it is unclear 
if the induced current from the atmosphere will have an 
appreciable effect on the interior dynamo, and its role 
sh ould be assessed with a deta iled, numerical model. As 
in (i Batygin fc Stevenson|20l"0 1 , we used a simple analyt- 
ical prescription tor the velocity field and the flow is con- 
flned to the radiative atmosphere so that that the RHS 
of the above equation vanishes in the (assumed) rigidly 
rotating interior. Nearly solid body rotation of the deep 
interior is predicted by the same arg ument that was used 
to reach this conclusion for Jupiter { ^Liu et al.|2008 ) . We 
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Fig. 2. — A series of representative pressure-temperature profiles 
of an evolved (t = 4.5Gyr) lMj„p planet, with T^ff = 1400K. 
The convective parts of the planet are plotted as thick red lines, 
while the radiative parts is plotted as thin purple lines. The four 
profiles correspond to e = 0%, 1%, 3% and 5% solutions. Note: the 
published ApJ version of this figure contains a plotting error. 
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calculation of u x Bdip, and implicitly in interpretation of 
GCMs of tidally locked planets. Once the electric scalar 
potential is known, the current, J, is determined from 
Ohm's law (including the induction cmf i.e. equation 3). 
Thus, the Ohmic dissipation per unit mass, P = J^/pcr 
everywhere in the planet can also be computed. 

As already discussed above, we can interpret the qual- 
itative characteristics of the velocity field from GCM 
results. Although the details of the flows that GCMs 
produce vary somewh at, depending on the assumptions 



implicit to the solver (Heng et al. 2011), most solutions 
produce a broad single jet in the atmosphere. The phys- 
ical mechanism for the jet's formation has been identi- 
fied as standing planetary-scal e Rossby and Kelvin waves 



(Showman and Polvani 2011). However, there exists a 



notable exception i n the literature, where ret rograde jets 



are also produced ( [Thrastarson fc Clio||2010 ). This may 
be of importance tor the Ohmic mechanism, since the 



importance 

atmospheric flows govern the induction of the interior 
current. To explore the effect of retrograde motions, we 
h ave computed the inter i or cur rents in a static model (as 



Batygin fc Stevenson] (20101) with three different pre 



scriptions tor the jets. iNamely, we considered, a single 
broad jet (v^ oc sin0), two counter- rotating jets, one in 
each hemisphere (w^ oc sin26'), as well as three jets, with 
retrograde equatorial flow (u^ oc sin 30). 

The results of these calculations are presented in side- 
view cross-sections of planets in Figure 1. In each of the 
three panels, the small arrows inside the planet show the 
induced current (actual model output). Additionally, the 
magnetic field lines as well as graphical representations 
of the currents and atmospheric flows are presented. Fig- 
ure lA shows the current induced by a si ngle jet. As in 
the results of jBatygin fc Stevenson (20101, loops are set 
up such that the radial current, which is induced at the 
equator flows towards the poles in the interior, and turns 
around in the atmosphere. If the current is induced by 
a double jet (Figure IB), its geometry changes consid- 
erably. Since the flows produce counter-acting electro- 
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Fig. 3. — A series of representative electrical conductivity pro- 
files of an evolved (t = 4.5Gyr) IMj^p planet, with T^ff = 1400-ft'. 
The ionization was derived from the the temperature-pressure pro- 
files, presented in Figure 2. The convective parts of the planet are 
plotted as thick red lines, while the radiative parts is plotted as 
thin purple lines. The four profiles correspond to e = 0%, 1%,3% 
and 5% solutions. Note: the published ApJ version of this figure 
contains a plotting error. 

motive forces, the interior current now flows from one 
hemisphere to the other. However, the hemispheric sym- 
metry is restored with three jets (Figure IC), and due to 
a particular choice of signs (namely, prograde flow at the 
equator in Figure lA and retrograde flow at the equa- 
tor in Figure IC), the direction of the interior current is 
reversed. 

Although, the interior current geometry changes con- 
siderably depending on the assumed flow, the result- 
ing heating rate does not. Indeed, a comparison of the 
Ohmic dissipation profiles in the three static models de- 
scribed above reveals that the ratio of the cumulative 
dissipation below the radiative-convective boundary to 
the total is the same to within ~ 10%. Consequently, 
in accord with most GCM results, we approximate the 
atmospheric flow as a single zonal jet in all following cal- 
culations. 

Our calculations of planetary interior structure and 
evolution employe d a descendant of t he Berkeley stel- 
lar evolution code ( Henyey et al.|1964 ). Over more than 
40 years of its lifetime, the model has been subject to 
a considerable amount of modification and improvement 
in its input physics and has been used extensively in the 
study of both extra-s olar and solar system giant planets 
'Pollack et al. 1996' ' ' ' 



Bodenheimer et al.| 2003 ; Batygin 
|etaLj|2009j) . A given evolutionary sequence starts at a 
radius of roughly 2Rjup and ends at an age of 5 Gyr. 

The program we used to compute the structural evo- 
lution of hot Jupiters assumes that the standard La- 
grangian equations of stellar structure apply. Energy 
transport is accomplished either by convection or radi- 
ation, as dictated by the Schwarzschild criterion. En- 
ergy sources within the planet include gravitational con- 
traction, cooling of the interior, and Oh mic dissipation. 
For the gaseous envelope, the interpolated [Saumon et al 
( jl995| ) equation of state is used. In the models with a 
core, the solid/liquid core has constant density. The im- 
perfections of the equation of state are expected to have 
a much smaller effect on the computed radii than the 
uncertainties in our model for Ohmic inflation. 
The atmosphere is taken to be gray. Pure molecular 
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Fig. 4. — A series of representative Ohmic heating profiles (en- 
ergy dissipation rate per unit mass) of an evolved (t = 4.5Gyr) 
lAfjup planet, with T^jf = 1400/C, corresponding to the e ^ 
temperature-pressure profiles showrn in Figure 2. The thick red 
lines denote the heating of the adiabatic interior, while the thin 
purple lines represent the heating of the radiative atmospheres. 
Note: the published ApJ version of this figure contains a plotting 



opacitie s are used in the radi ative outer layers of the 
planet (Freedman et al. 2008) while the high temper 



ature and pressure, opacities of [Alexander &: Ferguson 



by definition. However, significant drops in tempera- 
tures shortly above the photosphere can be o bserved in 
mo re sop histicated atmospheric models (e.g. [Spiegel et| 
al. (20091). Furthermore, if a given atmosphere exhibits 
aTxhermal inversion, the characteristic decrease in tem- 
perature at P ~ 10~^ bars provides a thin, electrically 
insulating layer. In other words, the null radial current 
boundary condition at the photosphere is likely to hold 
true for most atmospheres, whether or not an inversion 
exists. 

A series of representative Ohmic heating profiles (cor- 
responding to the e 7^ temperature-pressure profiles 
presented in Figure 2) are shown in Figure 4. Note 
that although the heating is maximum at the very sur- 
face, the layers that contribute to the expansion lie much 
deeper (i.e. below ~ 10 bars at an early epoch and be- 
low ~ 10^ — 10^ bars for evolved planets). The loca- 
tion of maximal heating corresponds to the atmospheric 
height where the electrical current turns around. Thus 
the higher up, electrically insulating layers of the atmo- 
sphere discussed in the previous paragraph, would expe- 
rience very little Ohmic heating. 

5. RESULTS: RADIAL EVOLUTION 

We computed the structural evolution of Ohmically 
heated hot Jupiters with the mass range spanning an 



|T994]) are used in the interior. In the atmosphere, it ^rder of magnitude between 0.23Mj„p and 3Mj„p and 



is assumed that dust grains contribute negligibly to the 
opacity. This may not be true but is again a smaller 
uncertainty than the other uncertainties in the m odel- 
ing. We refer the reader to Hubickyj et al. (2005) and 



Dodson-Robinson et al. ( 2008[) tor further reading. Oiven 
the approximate treatrnent of magnetic induction in the 
model, it is sensible to utilize a simple gray atmosphere. 
However, it is noteworthy that a different atmospheric 
model could modify our results on a quantitative level 
since the atmospheric boundary condi tion directly gov- 



erns the cooling rate of the interior (Guillot & Havel 



A series of representative temperature-pressure profiles 
of an evolved IMj^p planet are plotted in Figure 2. The 
four presented evolutionary sequences were started from 
the same initial condition, but were evolved with differ- 
ent amount of Ohmic heating, namely e = 0, e = 1%, 
e = 3% and e = 5%, as labeled. It is noteworthy that 
in planets with higher Tg/ / , Ohmic dissipation in the at- 
mosphere can lead to a formation of a small convective 
region, in the otherwise radiative atmosphere. Electri- 
cal conductivity profiles, corresponding to the presented 
temperature-pressure profiles are shown in Figure 3. 

At the planetary surface, defined as the Rosseland 
mean photosphere, the luminosity is composed of two 
parts: the internal luminosity generated by the planet, 
and the energy absorbed from the stellar radiation fiux 
and re-radiated (insolation). In our simulations, the at- 
mospheric circulation spans the region between 1 bars 



and the photosph ere, in accord with GCM results (Show- 
man et al. 2009). Additionally, in the heating calcula- 
tion, a boundary condition of zero radial electrical cur- 
rent is employed at the surface, as a result of signifi- 
cant decreas e in conductivity with heig ht in the upper 
atmosphere ( Batygin fc Stevenson||2010 ). This decrease 
in electrical conductivity is not exjalicitly present in our 
models, since the temperature at the photosphere is Tg// 



effective temperature range between lOOOiC and ISOOif . 
All but one of our models are core- less, and thus should 
be viewed as giving an upper bound on the radius that a 
planet of a given mass and Tf,f / may achieve. In our suite 
of simulations, we observed two distinct families of solu- 
tions that are characterized by either a monotonically 
decreasing or a monotonically increasing radius in time. 
A few representative evolution sequences are presented 
in Figure 5, with their masses and effective temperatures 
labeled. In the solutions where the radius is monoton- 
ically decreasing with time, the final answer (radius at 
t — 5Gyr) is largely independent of the initial condi- 
tion: the solution always asymptotically approaches the 
same "equilibrium radius." If the radius monotonically 
increases with time, however, the quantity of importance 
is the total integrated heating rate, which is intrinsically 
a function of the planetary radius. Consequently, these 
"unstable" solutions are unavoidably dependent on the 
initial conditions. 

Hot Jupiters can change their radii significantly pri- 
marily because of a change in the entropy of the deep 
interior, since the outermost several (radiative) scale 
heights constitute only a small fraction of the planet ra- 
dius. In the idealized limit where the heat escaping from 
the deep interior is small compared to Ohmic heating, 
the equation governing evolution for the deep interior 
reads 

|(i?grav + i?int)=7'- (19) 

If the planet is a polytrope and the internal energy is 
entirely thermal then the sum of internal and gravita- 
tional energies can be written as —kGM^/R wher e fc is a 
constant typica lly somewhat smaller than unity ( Chan- 
drasekhar 1957). In the case of highly degenerate bodies 
(e.g. Jupiter itself) the zero temperature part of the 
internal energy changes in such a way as to exactly can- 



Thermal Evolution of EGP Radii 



7 



R (Rjup) 




I Mj„pT,»-=l400K 



1 t (yr) 



10' I0» 10' 

Fig. 5. — A sample of evolution sequences of 0.5, 0.7 and lMj„p 
coreless planets at T^ff = 1400 and 1800K with t = 3%. The 
0.7 and IMj^p planets at 1400-K' planets asymptotically approach 
radii that are characterized by thermal equilibrium. The 0.5, 0.7 
and IMj^p planets at ISOOiC are all on unstable evolutionary 
paths which eventually lead to Roche-Lobe overflow. However, 
the IMj^p planet, however, is considerably more stable, owing to 
its more degenerate state. The evolution of the 0.5Mj„p planet 
at 1400-ff lies between the two regimes - its starting radius is the 
equilibrium radius. 

eel the change in gravitational energy as radius changes 



(Hubbard 1984) and the LHS above becomes only the 



time derivative of thermal energy to a good approxima- 
tion. (This invalidates the claim, often made but er- 
roneous, that the luminosity of degenerate planets is de- 
rived from contraction) . The bodies of interest to us have 
non-ideal thermodynamics and are not in the degenerate 
limit so no simple result holds. However, direct calcula- 
tion shows that it is still approximately true to replace 
the LHS above with something like d/dt{-kGA'P / R) 
with k of order unity, even though it is not correct to 
think of the energy as being derived solely from changes 
in gravitational energy. All of this discussion ignores 
possible effects arising from redistribution of the heavy 
elements. 

The transition from stable to unstable solutions can 
only be understood through detailed models, because it 
depends on the details of opacity and conductivity struc- 
ture, but the essential physics lies in the comparison of 
Ohmic heating at depth with the total radiative heat 
loss from the convective interior evaluated at the radia- 
tive/convective boundary. This heat loss scales as the 
adiabatic temparature gradient, which depends in turn 
on g, the gravitational acceleration. This is approxci- 
mately linear in mass because radius does not vary much, 
so Ohmic heating overwhelms heat loss from the interior 
once the mass is sufficiently small. In the limit where 
the Ohmic heating dominates the evolution of the deep 
interior, we expect that the characteristic timescale of 
inflation Ti„fi = R/{dR/dt) is (to order of magnitude) 
the ratio of GAP /R to Ohmic heating, or 



e Mj R leff 



(20) 



As can be seen in Figure 5, after a sufficient amount of 
time, unstable solutions enter a phase of runaway growth, 
with the instability driven by the fact that the cumula- 
tive heating rate is proportional to the surface area of 




lOOOK 



I400K 



Tefr(K) 



Fig. 6. — An R — T^ff diagram of evolved models of various 
masses (e = 1%) as well as the current aggregate of detected tran- 
siting hot Jupiters. The solid lines represent theoretical curves, 
corresponding to core-less models, while the dashed line depicts 
A/ = 0.5Af j„p models with a 2OM0 solid core. The observational 
data points are presented, along with their corresponding la error 
bars. The masses of the data are binned to fall between the the- 
oretical curves and are labeled by color and size as shown on the 
figure. 

the planet. This inevitably lead s to Roche-lobe ov erflow 
and evaporation of the planet (Laine et al. 2008). The 



timescale over which the planet can remain intact while 
on an unstable path is a function of the planetary mass 
and high-mass planets can, in principle, remain intact 
and grow slowly over many Gyr, while low-mass giant 
planets evaporate on a IGyr timescale. It is notewor- 
thy that the timescale will also likely be affected by the 
1/R slowdown of winds (see equation 16) and thus a pos- 
sible modification of the efficiency (which we kept con- 
stant) with a growing radius. Still, it is very likely that 
evaporation is unavoidable in certain circumstances. 

Low-mass planets can be stabilized against evapora- 
tion, i f they possess sufficiently massive high-met allicity 
cores ( |Bodenheimer et al.||2003 [Ibgui et"aL 2010[ ) . How- 
ever tne critical core mass is tar from being trivially 
small. In fact, a simulation of an O.SMjup planet with a 
p = 3g/cc, 2OM0, core at Tg// = 1800if reveals that it is 
only able to retain its envelope for 1 Gyr before over- 
flowing the Roche-lobe. Given that hot low- mass planets 



Papaloizou 


20071), while hot Jupiters 


by companions (|Lo Curto et al. 


2010 



cant fraction of hot, sub-giant planets, without compan- 
ions, may have been born asgiant planets and have since 
lost their gaseous envelope^ 

Within the context of our mechanism, the planetary 
radius is most clearly visualized as function of effective 
temperature rather than mass. Consequently, we choose 
to plot equal-efficiency, equal-mass lines on an i? — T^f / 
diagram rather than a conventional R — M diagram. 
Our results, along with the current aggregate of well- 
characterized transiting hot Jupiters, are illustrated in 
Figures 6-8. The evolved flnal model radii are plotted as 
points with spline interpolation between them and are la- 
beled according to their masses. The observational data[^ 

^ A similar idea has alre ady been proposed in the context of 
thermal atmospheric escape i jBaraflFe et al.|2005[ l. 

^ The data was acquired trom the Kxtrasolar Encyclopedia 
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Batygin Stevenson Bodenheimer 




Fig. 7. — Same as Figure 2, but with e = 3% 

are plotted as points with their corresponding la error 
bars. The data are binned in mass, to faU between our 
model masses, and are labeled by color. As a result, if 
the Ohmic mechanism is indeed correct, and unmodeled 
effects such as var iable opacity a nd age make a small 



contribution (Burrows et al 
a given color "should tali be 
same color. 



2007 1, all observed data of 
ow the model curve of the 



6. DISCUSSION 

As is clearly demonstrated, if the approximations in- 
herent to our model are valid, Ohmic dissipation, along 
with the likely presence of high-metalicity cores, is able 
to account for all of the currently observed radius anoma- 
lies. Furthermore, the self-similar family of model curves 
reproduces the characteristic increase and overturn of the 
radii in the effective temperature region between 1200X 
and 1800/f . This feature is unique to the Ohmic dis- 
sipation mechanism and is a consequence of ionization. 
Recall that the primary source of electrons in the atmo- 
sphere is alkali metals, with K playing a dominant role. 
Under densities that are characteristic of hot Jupiter 
atmospheres, K has an ionization temperature around 
~ 1800-?^. This explains the rapid increase in Ohmic 
heating with rising Tg// around ~ 1400K. This increase 
comes about through a Boltzmann factor in the Saha 
equation, and will therefore have a roughly exponential 
dependence. We stress that this is uniquely characteris- 
tic of the Ohmic mechanism, since other proposed mech- 
anisms will exhibit an algebraic dependence on T^jf (see 
Laughlin et al. (2011 ) for a statistical comparison of var- 
ious inflation mechanisms with data). 

As the temperature nears ^ 1800K, the ionization lev- 
els begin to saturate and the radiative/convective bound- 
ary moves deeper, so the radii don't grow as quickly with 
Tf,ff. However ionization of Na begins to contribute 
strongly above 2000X. This renders all but the most 
massive of our models unstable. This is suggestive, given 
that out of a sample of 87 transiting Hot Jupiters, only 
four have been detected with Te// > 2100isr, although 
the mechanism by which Hot Jupiters halt their migra- 
tion at small radii is still not fully understood. 

It is noteworthy that our low-mass models signifi- 
cantly overestimate the data, to the extent where small 
(~ 2OM0) cores are insufficient to explain the observed 

( jhttp:/ / www.exoplanet .eu| ) 



Fig. 8. — Same as Figure 2, but with e = 5% 

radii. This is not particularly surprising, given that the 
inflation is somewhat enhanced by a shallower adiabatic 
temperature gradient, g/cp, allowing for a larger frac- 
tion of the generated heat to be deposited into the con- 
vective interior, and a diminished gravitational energy, 
that must be overcome by Ohmic heating. The discrep- 
ancy between the model and the data is either indicative 
of massive cores, or more likely curbed efficiency of the 
Ohmic mechanism at smaller mass. Both are certainly 
plausible, and the latter can be a consequence of dimin- 
ished strength of the magnetic fi eld, which is di ctated by 



the intrinsic heat-flux | Christensen et al.||2009| |. 

As the parameter space available to hot Jupiters is, to 
first order, delineated, the models presented in this work 
make two predictions. First and foremost, the radii of 
all hot Jupiters should fall below the theoretical curves 
presented in Figures 6-8. Second, a fraction of the hot 
sub-giant planet population may have originated from 
the giant planet population. Since this process can take 
place many Gyr after the birth of the star, Ohmically 
evaporating planets should be detected, given a suffi- 
ciantly large data set. These predictions will become 
testable, as the known aggregate of transiting exoplan- 
ets continues to grow. Accordingly, additional data will 
allow us to further understand the diversity of planetary 
bodies and to better define the Solar System's place in 
our galactic neighborhood. 

We would like to conclude with a discussion of future 
prospects for improvement of the model. Let us begin 
by recalling the simplifications inherent to the approach 
presented here. In the structural calculation, we used 
a gray, one- dimensional atmosphere. Obviously this is 
not intended to be a quantitatively good representation 
of real atmospheres and improved treatment of the tem- 
perature (and thus conductivity) profile, accounting for 
asymmetry in insolation between the day and the night 
side as well as heat transport by global circulation is de- 
sired in the future. Indeed the considerable difference in 
interior cooling rates sh own here, and those presented by 



Guillot fc Havel 
ing the atmosp 



(2011|) suggests that the effect of chang- 
oundary condition on structural 
That said, the inability 



)neric 

evolution may be considerable, 
of current state-of-the-art atmospheric models to match 
observations in a satisfactory manner suggests that much 
additional effort must be applied towards improvement 
of our understanding of extra-solar planet atmospheres 
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for implementation of a complex atmosphere into our 
model to be meaningful. In the heating calculation, our 
analytical parameterization of the circulation should be 
replaced by a self-consistent solution of the global cir- 
culation that fully accounts for magneto-hydrodynamic 
effects. The induced current in the interior would be a 
natural result of such a calculation. At the same time, 
the effect of the induced current on the interior dynamo 
should be assessed. Although in this study we assumed 
that the interior dynamo is largely unaffected by the cur- 
rents induced in the atmosphere, in reality this may not 
be the case, since the induced current can in principle be 
comparable to the characteristic current of the interior 
dynamo. However, given the enormous computational 
cost of such an endeavor, it is not surprising that a cal- 
culation of this sort has not yet been done. 

In general, future efforts in improvement of the 
Ohmic heating model should rely on both, modeling 
and observational constraints. After all, observations 
of primary and secondary eclipses can often yield 
information about th e atmospheric temperature profile 
Knutson et al. '2008') as well as elemental abundances 
Line et al.,,2U10j . At the same time, ground-based 



high-dispersion sp ectroscopy can yiel d an estimate of 
the wind velocity (Snellen et al. 2010). In certain cases 
(e.g. WASP-12b), the planet's proximity to its host 
star results in considerable interactions between the 
planetary magnetosphere and the stellar wind. This 
can lead to formation of a bow-shock, which ultimately 
produces an observable signature, that manifests itself 
as an early ingress of the transit in near-UV frequency 
range. Quantitative modeling of this process may lead to 
meaningful constraints on the planetary magnetic field 
(Vidotto et al. 2010). As a result, an observationally 
guided improvement of the Ohmic inflation model will 
surely prove to be a useful venture, as quantitative 
enhancement of the first-order results presented in this 
paper will allow for a rigorous statistical comparison 
between the model and an expanded data set. 
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